The extensive age gradient of the Carina dwarf galaxy 
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ABSTRACT 

The evolution of small systems such as dwarf spheroidal galaxies (dSph) is 
likely to have been a balance between external environmental effects and inter- 
nal processes within their own relatively shallow potential wells. Assessing how 
strong such environmental interactions may have been is therefore an important 
element in understanding the baryonic evolution of dSphs and their derived dark 
matter distribution. 

Here we present results from a wide-area CTIO/MOSAIC II photometric 
survey of the Carina dSph, reaching down to about two magnitudes below the 
oldest main sequence turn-off (MSTO). This data-set enables us to trace the 
structure of Carina in detail out to very large distances from its center, and as a 
function of stellar age. 

We observe the presence of an extended structure made up primarily of an- 
cient MSTO stars, at distances between 25'-60' from Carina's center, confirming 
results in the literature that Carina extends well beyond its nominal tidal radius. 

The large number statistics of our survey reveals features such as isophote 
twists and tails that had gone undetected in other previous shallower surveys. 
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This is the first time that such unambiguous signs of tidal disruption have been 
found in a Milky Way "classical" dwarf other than Sagittarius. 

We also demonstrate the presence of a negative age gradient in Carina directly 
from its MSTOs, and trace it out to very large distances from the galaxy center. 
The signs of interaction with the Milky Way make it unclear whether the age 
gradient was already in place before Carina underwent tidal disruption. 

Subject headings: galaxies: dwarf — galaxies: evolution — galaxies: individ- 
ual(Carina) — galaxies: stellar content — galaxies: structure — Local Group 



1. Introduction 

The dwarf spheroidal galaxies (dSph) surr ounding the Milky W ay (MW) are the galaxies 
studied in better detail than any other (see iTolstoy et al.l 120091 for a recent review), yet 
many fundamental questions about their nature remain unanswered. Are they the most 
dark matter dominated galaxies we know of, as their large apparent dynamical mass-to-light 
ratios would suggest, or is tidal disruption a contributing factor ? What drives the observed 
variety of star formation and chemical enrichment histories ? Were dSphs born as we see 
them today or are they the descendants of larger and intrinsically different systems? 

The evolution of small systems such as dwarf galaxies is likely to be particularly sus- 
pectible to the environment that they have experienced as well as being influenced by their 
own low potential well. However, the relative importance of both these factors has yet to be 
established. 

Among the MW satellites, an example of a dSph whose evolution may have been strongly 
affected by the environment is given by Carina. This galaxy is the one "classical" dSph where 
most hints of tidal disruption have been found, such as the presence of probable members 
well beyond the nominal tidal radius, a rising velocity dispersion profile and a velocity 

and references 



2006 



shear reminiscent of tidally disrupted dwarf galaxies (see iMunoz et al.l 
therein). However, no report of such kinematical features is given in iKoch et al.l ( 120061 ) and 
Walker et al.l ( 120091 ). so that the hypothesis of tidal disruption awaits confirmation. 



Carina is a particularly interesting object among the Local Group dwarf galaxies also 
because it is the only clear example of a galaxy with well-separated star formation episodes, 



Bono et al. 


2010). The effects of its i 


as well ( 


Koch et al. 


2008 




Venn et al. 



Lemasle et al.ll2012l ). It is unclear at present whether the influence of the MW may have 
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induced the main SF episodes (e.g. see simulations by lPasetto et al.ll201l[ ) as episodic SFHs 
can be qualitative ly explained also by models of low mass dSphs treated in isolation (e.g. 
Revaz et al-lboogh . 



Another feature that can provide information on the evolution of a small galaxy is 
the spatial distribution and kinematics of its stellar population mix. Negative age and 
metallicity gradients are found in several of the Local Group dSphs, both around the M W 
JHarbeck et aliboOli llblstov et al.lbooi IBattadia et al.lbood . boill ; Ide Boer et d]l2012aH ) 
and in isolation (e.g. Bernard et al .12008 ); the latter finding suggests that internal effects 
are the main drivers in setting up the observed gradients. In Carina, indications of an age 
gradient have been found from th e spatial distribution of red clump (RC) stars as compared 
to horizontal branch (HB) stars (jHarbeck et al.ll200ll ). Although this is similar to the be- 
haviour of other MW dSphs with multiple stellar components, for Carina there have been no 
findi ngs of different kinematics am ong the stars belonging to different age/metallicity ranges 
(e.g. IWalker fc Penarrubial I2OIII). unlike what observed in Fornax, Sculpt or and Sextans 
jTolstov et alil2004i IBattadia et al.lbood boill IWalker fc Penarrubialhoilh . Interestingly, 
it has been shown that strong interaction with the MW is a viable mechanis m for erasing 
differences in the kinematics of multiple stellar components ( jSales et al.ll2010l ). 



Several of the observed properties of Carina could therefore be explained in a framework 
of strong interaction with the MW. It is then particularly interesting to fully characterize 
the stellar population of this object and to put on a secure ground the possible existence of 
tidal features. 

Here we present results from a deep, wide-area photometric survey of Carina, which 
yields major improvements with respect to previous studies of the morphology and spatial 
distribution of Carina's stellar population mix. 



2. Observations 



Deep optical photometry covering the Carina dSph was obtained using the CTIO 4-m 
MOSAIC II camera over 6 nights in December 2007, as part of observing proposal 2007B- 
0232 (PI: M. Mateo). Observations in the B and V bands were obtained in 8 pointings 
around the centre of Carina, covering a total area of ~2 deg 2 . The observing strategy 
was to obtain multiple long (^600s) exposures for each pointing, which have been stacked 
together to obtain deep photometry. In addition, short (30s) exposures were also obtained, 
to sample bright stars that are saturated in the deep images. To en sure accurate photometric 
calibration of the data set, observations of Landolt standard fields (jLandoltlll992[ ) were made, 
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covering a range of airmass and colour. 



First stage image processing steps included standard bias and overscan-correction, trim- 
ming to the reliable active detector area, and flatfielding via creation of master flats from a 
dithered set of B-band and V-band twilight sky exposures. The flatfielding step also corrects 
for internal gain variations between the detectors. 

Prior to deep stacking of the individual exposures, detector-level catalogues were gener- 
ated for each science image to refine the astrometric calibration and also to assess the data 
quali ty. For astrometric calibration, a Zenithal polynomial projection (IGreisen fc Calabretta 
20021 ) was used to define the World Coordinate System. A third-order polynomial includes 
all the significant telescope radial field distortions leaving just a six-parameter linear model 
per detector to completely define the astr ometric transform ations. The TwoMicronAll Sky 
Survey (2MASS) point-source catalogue (jCutri et al.l 120031 ) was used for the astrometric 
reference system. 

The common background regions in the overlap area from each image in the stack were 
used to compensate for sky variations during the exposure sequence and the final stack 
included seeing weighting, confidence (i.e. variance) map weighting and clipping of cosmic 
rays. 

Catalogues for each de ep stacked image were then generated using standard aperture 
photometry t echniques (e.g. Ilrwin et al.ll2004l ) and cross-matching with published photome- 
try of Carina (IStetson 1120001 ) was used to directly photometrically calibrate the central fields. 
Outlying fields were calibrated with respect to the central region using the overlaps between 
adjacent pointings. 



3. Results 



3.1. Age Gradient 



While the exact age of Carina's SF episodes varies among the various works in the 
literature, there is consensus about Carina experiencing at least three main SF episodes: an 
"an cient" one, around 11-13 Gyr ago; a main episode at intermediate ages (around 2.5-4 Gyr 



ago iBono et al.l |2010| . but ~7 Gyr ago according to Hurley-K eller et al. 1998), af ter a long 



quiescent period; a relatively recent episode, >1 Gyr ago (e.g. iMonelli et al.ll20031 ). 



The way these "ancient" , "intermediate-age" and "young" stellar components are spa- 
tially distributed has not been studied in detail yet. Notably, notwithstanding its complex 
SFH, Carina has a very narrow red giant branch (RGB). This makes photometry of RGB 
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stars unsuitable for exploring the existence o f age gradients in t his galaxy, even though such 



spatially extended data-sets exist f~1 0deg 2 IMunoz et al 



2006). Indications for a negative 



age gradient in Carina were found by lHarbeck et al. I feooih over a relatively small radial 
range (30'), using photometry down to ~lmag below the HB. The shallow photometry al- 
lowed only the use of indirect age indicators for tracing the gradient, i.e. the RC for the 
young/intermediate-age stars and the HB for the ancient stars. 

Figure [T] shows the Colour-Magnitude Diagram (CMD) resulting from our photometric 
survey, split over a selection of elliptical radii R. This spatially extended photometric data- 
set reaches down to ~2 magnitudes below the oldest MSTO, allowing us to trace possible 
age gradients directly from the presence or absence of the MSTOs belonging to the various 
SF episodes and over a large area. 

In the innermost region (R < 10') all the main features of the Carina stellar population 
are present: the thin RGB, the RC, the red and blue HB, the ancient MSTO and sub-giant 
(SGB), the intermediate-age MSTO and SGB, and younger MS stars at B — V~0 and 21< V 
<_23 (for a detailed discussion on the various features of the CMD see e.g. iHur ley- Keller et al. 



19981 ). At larger radii (10' <R< 25') the young stars are barely visible, while the features 
containing intermediate-age and ancient stars are still well populated. In the outer parts 
(R> 25'), the MSTO of the intermediate-age population disappears, while the ancient stars 
are still visible, as witnessed by the presence of the old MSTO and SGB. 

The above direct detectio n of an age gradient in Carina, as traced from MSTO stars, 
confirms the early findings of lHarbeck et al.l (120011 ) and extends them over a much larger 
radial range. 

We note that the presence of a sizeable population of ancient MSTO stars in the last 
two panels of Figure [T] (R> 25' ) shows that Carina e xtends well beyond its nominal King 
tidal radius (r t =28.8', llrwin fc Hatzidimitrioul [1995| , hereafter IH95). This puts on a se- 
cure ground the conclusions by IMunoz et al.l (120061 ) based on the spatial distribution and 
kinematics of a smaller sample RGB stars. 

Both the presence of the age gradient and the large extent of Carina are also illustrated 
by Figure |2j where we compare to a control field the luminosity function (LF) in bins of R 
of stars in the region of the MS, MSTO and SGB. Owing to the large extent of Carina it 
is difficult to accurately identify a region containing only contaminants (i.e. MW stars and 
unresolved galaxies) to use as a control field. We chose to use a sub-region of the outermost 
pointing on the North-East (see Figure |3]) because it has the lowest density of objects among 
the various pointings; the significantly different shape of the control field LF from the two 
outer elliptical zones confirms that it mostly contains contaminants and demonstrates the 



-6 - 



presence of appreciable numbers of Carina members out to 1 degree radius. Furthermore, 
the similarity of the LFs of the two outer zones and the absence of obvious intermediate-age 
and young stars beyond ~25' strongly suggest this extended structure is made up primarily 
of ancient stars. 



3.2. Extended Spatial Structure 

Figure [3^ illustrates the spatial distribution of the overall Carina stellar population after 
applying magnitude cuts of B<25.1, V<25.5 - corresponding to ~50% completeness levels 
- to counter the variable depth of the pointings due to differences in total exposure time. 
To minimize the presence of contaminants, we isolate the main features of Carina's overall 
stellar population on the CMD via selection in magnitude and color. Figure [3b,c,d show 
how the spatial distribution varies with age. 

The isodensity contours in each panel were derived from maps of stellar counts (Hess 
diagrams) on a 2' x 2' pixel grid, from which we subtracted the density of contaminants, 
^cont- This was calculated as the median of the pixel values in the control field (shown in 
Fig. [3] and discussed in 13.11) . The iso-density contours are then plotted in units of noise 
(a), which we determine as the r.m.s. of the pixel values in the contaminant-substracted 
control field. Should the derived S cont be overestimated because of a leakage of Carina's 
stars, the main effect will cause Carina to appear smaller than it is and to underestimate 
the significance of the contours; hence the extent of Carina found here shold be considered 
as a "lower limit" and the significance of the outer contours as conservative. 

Figure Eh shows that the main body of Carina appears of regular morpology, although 
with slightly boxy iso-density contours, out to i?~18-20' ~ 0.6-0.7 r t . Beyond this radius, the 
contours become elongated and distort, with very extended tails emerging from both sides 
of the main body and a well-defined over-density of stars at (£,?7)~(— 60',— 30'). The much 
larger number statistics of this study with respect to studies in the literature makes it possible 
to observe isophote elongation and twists; these had not be detected in previous, shallower 
photometric su rveys of Carina, eye n though optimized to minimize the contamination from 



MW stars (e.g. Munoz et all 120061 ) . 



Fitting of ellipses as function of R to the map of stellar counts confirms the change in 
shape and orientation of the contours, yielding a P.A. and ellipticity e that are constant 
around a P.A. = 53.5° ±1.0° and e = 0.29±0.01 at R < 18' , and then increase, yielding a 
P.A.= 62.5° ±0.5° and e = 0.35±0.02 at R > 18' ; the latter values tend to the literature 



values of PA.=65° ±5° , e =0.33±0.05 fllrwin fc Hatzidimitrioulll995f ) 
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We note also the presence of a break in the surface number count profile of Carina's stars 
lying along the projected major axis of the system (Figure H]). The break occurs in both N-E 
and S-W halves of the galaxy at around 20'from the center, where the iso- density contours 



start becoming elonga t ed, and at a similar r adial distance to that found by IMajewski et al. 
J2000h : IWalcher et all fcpQ3h : iMunoz et all tood e.g.) in shallower samples. 

The youngest Carina stars do not extend as far out as the intermediate-age and ancient 
ones, but signs of isophote elongation are arguably also visible in the spatial distribution 
of these stars (Figure [3] b). This would set the last pericentric passage less than 1-2 Gyr 
ago. The intermediate-age and ancient population both show distorted outer isophotes 
(Figure [3] c,d). Notably, as visible in the figure and confirmed by ellipse fitting of the Hess 
diagrams, the shape parameters of the contours are subtly different for these two stellar 
populations of different ages and may impact on models trying to reconstruct the orbital 
history of Carina. 



Discussion and conclusions 



We have used CTIO/MOSAIC II photometry to study the morphology and spatial 
distribution of Carina's stellar population mix. 

The combination of wide-area and photometric depth, reaching down to ~2 mag below 
the oldest MSTO, allows for a direct detection of an extensive age gradient in Carina, traced 
from the presence/absence of the MSTOs belonging to the various SF episodes. 

We find a clear presence of ancient MSTO stars at R > 25' , which demonstrates that 
Carina extends well beyond its literature nominal tidal radius, confirming results from pre- 
vious works. The increased number statistics from this study reveals the existence of faint 
features such as isophote elongation and twists, that could not be detected in previous shal- 
lower photometric surveys. Such features are even visible when separating the stars in age. 
We confirm the presence of a clear break in the surface brightness profile at radii >20 ' and 
observe an over-density of stars in the S-W half of Carina, at ~70' from the center. 

Features such as breaks in the surface brightness profile, together with variations in the 
elongation and direction of the isophotes, have been detected in o ther objects that are good 
candidates for ongoing tidal distruption such as NGC205 and M32 (IChoi et al.ll2002f). and are 



predicted by N-bqdy simulations of tidally disrupted dwarf galaxies (e.g. 



Choi et al. 2002 



Munoz et al.l 120081 ; iPenarrubia et al.l 120091 ) . These morphological traits strongly suggests 



that Carina is being tidally disrupted. 
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At a distance of ~100kpc the low Galactocentric line-of-sight velocity of Carina (~7 
km/s) places it either close to apocenter or pericenter. Given the signs of tidal disturbance 
and the current best estimate of the proper motion (IPiatek et al.l 120031 ) . Carina is most 
probably close to apocenter with an orbital period of between 1-2 Gyr and hence the last 
pericentric passage likely occurred <1 Gyr ago. The presence of si gns of tidal disruption a t 



apocenter is consistent with the timescales from the simulations of iPenarrubia et al.l ( 120091 ). 
where tidal tails can remain clearly visible at projected radii larger than the break radius 
until up to ~15 dynamical crossing times, i.e. between 1-2 Gyr in this case. 

The existence of signs of tidal disruption complicates the interpretation of the observed 
age gradient, as it is not obvious whether it was (in part) already in place before Carina was 
tidally disturbed. Related questions are whether Carina showed distinct kinematic signatures 
for the different populations prior to tidal interactions and how these may have been affected 
by close pericentric passages. 

This also impacts on the derived dark matter content and distribution. Treated as an iso- 
lated equil ibrium system, Car ina has a large dynamical mass-to- light ratio (e.g. M/Ly ~100 



(M/L)v,0, [Walker et al.l 120071 ). as do most other "cl assical" dSph s. An attempt to model 



Carina as a tidally perturbed system was made by IMunoz et al.l ( 120081 ). who found that 
the surface brightness and internal kinematic properties of Carina are compatible with a 
mass-follows-light model still containing large amounts of dark matter (M/L~40). The 2D 
distribution resulting from such model is remarkably regular, in contrast with what is ob- 
served he re over the s a me sp atial region. However, it should be considered that the best 
model of IMunoz et al.l ( 120061 ) was aiming at reproducing the spatial distribution from the 
data available at the time, too shallow to reveal the signs of tidal distorsion seen in our 
survey. 

Knowledge of the distorted structure of the galaxy from this deep survey, combined 



Munoz et al. 


2006; 


Walker et al. 


2009) 



(e.g. 



Koch et al. 2006 



the impact of tidal disruption on the derived mass-to-light ratio - and dark matter content. 

Tidal disruption of dwarf galaxies has been detected in some M 31 satellites (e.g. 
Choi et all [2002i iGeha et alJbood : iMcConnachie k Irwiafeood ). However, around the MW, 
apart from the clear case of Sagittarius, the quest has proven difficult due to the combi- 
nation of the large extent on the sky of MW dSphs and their low surface brightness. The 
question therefore remains if tidal disruption is common among the MW dSphs and ultra- 
faint dwarfs, or wheth er Carina belongs to a "minority" of tidally disrupted objects such as 
the Sagittarius dwarf (llbata et al.lll994l ). Hints of tidal disruption have been found in some 



of these other systems, e.g. Ursa Minor (e.g. iMartmez-Delgado et al.l l200ll ; iPalma et al. 
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20031 ). Leo I flMateo et al.ll2008f ). Ursa Major II (IMunoz et alJboiOf l. In contrast, Draco, the 



only other dSph surveyed with photomet ric data as deep an d wide those used here, 

shows a "flawless", regular morphology (jSegall et al.l 120071 ) and an enourmous M/L ~300 
falevna et aDl2002h . 



The research leading to these results has received funding from the European Union 
Seventh Framework Program (FP7/2007-2013) under grant agreement number PIEF-GA- 
2010-274151. The authors acknowledge the International Space Science Institute at Bern for 
their funding of the team "Defining the full life-cycle of dwarf galaxy evolution: the Local 
Universe as a template" . 

Facilities: Blanco (Mosaic II) 



REFERENCES 

Battaglia, G., Tolstoy, E., Helmi, A., et al. 2006, A&A, 459, 423 
Battaglia, G., Tolstoy, E., Helmi, A., et al. 2011, MNRAS, 411, 1013 
Bernard, E. J., Gallart, C, Monelli, M., et al. 2008, ApJ, 678, L21 
de Boer, T. J. L., Tolstoy, E., Hill, V., et al. 2012, A&A, 539, A103 
de Boer, T. J. L., Tolstoy, E., Hill, V., et al. 2012, A&A, 544, A73 
Bono, G., Stetson, P. B., Walker, A. R., et al. 2010, PASP, 122, 651 
Choi, P. I., Guhathakurta, P., & Johnston, K. V. 2002, AJ, 124, 310 

Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, "The IRSA 2MASS All-Sky Point 
Source Catalog, NASA/IPAC Infrared Science Archive. 

Geha, M., Guhathakurta, P., Rich, R. M., & Cooper, M. C. 2006, AJ, 131, 332 

Greisen, E. W., & Calabretta, M. R. 2002, A&A, 395, 1061 

Harbeck, D., Grebel, E. K., Holtzman, J., et al. 2001, AJ, 122, 3092 

Hurley-Keller, D., Mateo, M., & Nemec, J. 1998, AJ, 115, 1840 

Kleyna, J., Wilkinson, M. I., Evans, N. W., Gilmore, G., & Frayn, C. 2002, MNRAS, 330, 
792 



Koch, A., Grebel, E. K., Wyse, R. F. G., et al. 2006, A J, 131, 895 
Koch, A., Grebel, E. K., Gilmore, G. F., et al. 2008, AJ, 135 
Ibata, R. A., Gilmore, G., & Irwin, M. J. 1994, Nature, 370, 194 
Irwin, M., & Hatzidimitriou, D. 1995, MNRAS, 277, 1354 
Irwin, M. J., Lewis, J., Hodgkin, S., et al. 2004, Proc. SPIE, 5493, 411 
Landolt, A. U. 1992, AJ, 104, 340 

Lemasle, B., Hill, V., Tolstoy, E., et al. 2012, A&A, 538, A100 

Majewski, S. R., Ostheimer, J. C, Patterson, R. J., et al. 2000, AJ, 119, 760 

Martmez-Delgado, D., Alonso- Garcia, J., Aparicio, A., & Gomez-Flechoso, M. A. 2001, ApJ, 
549, L63 

Mateo, M. L. 1998, ARA&A, 36, 435 

Mateo, M., Olszewski, E. W., & Walker, M. G. 2008, ApJ, 675, 201 

McConnachie, A. W., & Irwin, M. J. 2006, MNRAS, 365, 1263 

Monelli, M., Pulone, L., Corsi, C. E., et al. 2003, AJ, 126, 218 

Munoz, R. R., Majewski, S. R., Zaggia, S., et al. 2006, ApJ, 649, 201 

Munoz, R. R., Majewski, S. R., & Johnston, K. V. 2008, ApJ, 679, 346 

Munoz, R. R., Geha, M., & Willman, B. 2010, AJ, 140, 138 

Palma, C, Majewski, S. R., Siegel, M. H., et al. 2003, AJ, 125, 1352 

Pasetto, S., Grebel, E. K., Berczik, P., Chiosi, C, & Spurzem, R. 2011, A&A, 525, A99 

Piatek, S., Pryor, C, Olszewski, E. W., et al. 2003, AJ, 126, 2346 

Peharrubia, J., Navarro, J. F., McConnachie, A. W., & Martin, N. F. 2009, ApJ, 698, 222 
Revaz, Y., Jablonka, P., Sawala, T., et al. 2009, A&A, 501, 189 

Sales, L. V., Helmi, A., & Battaglia, G. 2010, Advances in Astronomy, 2010, article id. 
194345 



- 11 - 



Segall, M., Ibata, R. A., Irwin, M. J., Martin, N. F., & Chapman, S. 2007, MNRAS, 375, 
831 

Smecker-Hane, T. A., Stetson, P. B., Hesser, J. E., & Vandenberg, D. A. 1996, From Stars 
to Galaxies: the Impact of Stellar Physics on Galaxy Evolution, 98, 328 

Stetson, P. 2000, PASP, 112, 925 

Tolstoy, E., Irwin, M. J., Helmi, A., et al. 2004, ApJ, 617, L119 
Venn, K. A., Shetrone, M. D., Irwin, M. J., et al. 2012, ApJ, 751, 102 
Walcher, C. J., Fried, J. W., Burkert, A., & Klessen, R. S. 2003, A&A, 406, 847 
Walker, M. G., Mateo, M., Olszewski, E. W., et al. 2007, ApJ, 667, L53 
Walker, M. G., Mateo, M., & Olszewski, E. W. 2009, AJ, 137, 3100 
Walker, M. G., & Peharrubia, J. 2011, ApJ, 742, 20 
Tolstoy, E., Hill, V., & Tosi, M. 2009, ARA&A, 47, 371 



This preprint was prepared with the A AS IATgX macros v5.2. 



-12- 



' — 1 — 1 — i — 1 — 1 — 1 — 1 — r 
R [arcmin] < 10 



T 1 1 1 1 1 1 1 1— , 1 ,— , 1- 

10 < R [arcmin] < 25 / ' • 



o 

> 



18 
20 
22 

24 

26 
18 

20 

22 

24 
26 



H 1 1 1 1 1 1 1 1 1 1 1 1 1 h 



H 1 H 1 1 1 1 1 1 1 1 1 1 1 1- 



25 < R [arcmin] < 45'- 



R [arcmin] > 45 



> 



-0.5 0.0 0.5 
(B-V)jc 



1.0 -0.5 



0.0 0.5 
(B"V) JC 



1 .0 



Fig. 1. — Colour-magnitude diagram of stellar objects along the line-of-sight to the Carina 
dSph, over different ranges of projected elliptical radius (see labels). The elliptical annuli 
have P. A. and e values from IH95. All the panels contain the same number of stars (15000) 
to enable a comparison of the relative strength of the features in the CMD. The boxes in 
the top- left panel show the cuts in magnitude and color applied to select young (blue), 
intermediate-age (yellow) and old (red) stars for Figure [3] 
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Fig. 2. — Luminosity function of Carina's stars, normalized to the peak value and plotted 
in bins of elliptical radii (see legend). The stars are selected to have 21 <V<25.5, B-V<0.6. 
The control field corresponds to the region with rectangular coordinates 60' < £ < 80'and 
20' < 7] < 40' (see Fig. E}. 
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Fig. 3. — Spatial distribution of stars in Carina for the overall stellar population, young stars, 
intermediate age population, old stars (see labels and Figured]). The overlaid iso-density 
contours were derived from the contaminant- density subtracted Hess diagram obtained on 
a 2'x2'pixel grid, and show levels of 1,2, 3,4, 5, 6, 7, 8, 9, 10a, 15, 20a... (where a is the noise, 
see text for details). The values of a are 2, 0.51, 0.24, 0.24 stars/pixel for panels a,b,c,d, 
respectively, while the values of the subtracted S cont are 3.08, 0.007, 0.01, 0.01 stars/pixel; 
the region used for determining a and S cont is enclosed by the dotted lines. The ellipse 
indicates the nominal tidal radius (parameters from IH95; center coordinates from Mateo 
1998). The horizontal dashed lines are meant to guide the eye when comparing the direction 
of the isophote elongation for the various populations. 
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Fig. 4. — Surface number count profile of Carina, derived from the stars in Figure Et) that 
lie on a 20' wide slice running along a P.A.= 65°. Filled and open squares refer to the North- 
East and South-West part of the galaxy, respectively. The contaminant density subtracted 
from the profile is indicated by the dotted line; the arrows show the radial range occupied 
by the region for calculating S cont . We note the change of slope at R~20' present in both 
halves of the galaxy, and the overdensity at R~70' in the S-W half, corresponding to the 
overdensity visible in Figure Ek) at (£,7/)~(— 60',— 30'). 



